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Abstract 


Radar  scattering  by  an  array  of  surface  land  mines  is 
studied.  The  array  Is  considered  to  be  a  random  perturbation  of 
a  uniform  array.  The  analytical  evaluation  of  the  expectation  Is 
performed  for  this  problem  under  specified,  but  realistic,  assump¬ 
tions.  The  analytic  expression  for  the  expectation  contains 
coherent  and  Incoherent  array  factors,  each  of  which  are  summations 
of  terms  that  are  weighted  by  the  effect  of  the  antenna  pattern. 

By  curve  fitting  the  antenna  gain  pattern  with  exponential  functions 
In  elevation  and  azimuth  the  resultant  series  are  summable.  The 
resultant  expression  for  normalized  signal  to  clutter  ratio  displays 
the  coherent  contributions  from  the  radar  system  parameters  (such 
as,  beamwidths,  frequencies,  depression  angle  and  pulse  width)  In 
conjunction  with  the  randomness  of  mine  placement  and  the  clutter 
distribution. 


TABLE  OF  CONTENTS 


Section 

ABSTRACT 

I  INTRODUCTION 

II  RESULTS  AND  OBSERVATIONS 

III  METHOD  OF  ANALYSIS 

A.  RADAR  SCAHERING  FUNDAMENTALS 

B.  ARRAY  FACTOR  DERIVATION 

C.  ARRAY  FACTOR  SUftlATIONS 

D.  SIGNAL  TO  CLUTTER  RATIO 

E.  GENERAL  CONCLUSIONS 

APPENDIX  I  SUMMATION  FORMULAE 

APPENDIX  II  SCAHERING  FACTORS  FOR  SMALL  MINE  FIELDS 

APPENDIX  III  COMPUTER  PROGRAMS 


REFERENCES 


I. 


INTRODUCTION 


Our  goal  In  this  work  Is  to  Investigate  the  process  of 
detection  of  randomly  spaced  mines  as  a  function  of  the  deterministic 
and  statistical  parameters  of  the  radar,  mine  distribution  and 
ground  clutter.  The  main  reason  for  such  an  undertaking  is  to 
provide  information  for  optimal  design  of  radar  systems  and  their 
implementation  in  the  mine  field  mapping  process. 

As  is  common  for  surface  mapping  radars  the  clutter  echoes 
are  dominant  over  system  generated  noise  and  a  key  parameter  in 
defining  the  radar  system  performance  is  the  echo  signal  power  to 
clutter  ration  S/C  [1;  Chap.  131  To  evaluate  this  ratio,  as 
indicated  in  the  proposal,  requires  the  computation  of  the  time- 
average  radar  echo  power  as  contributed  by  the  randomly  placed 
scatterers  in  the  range-azimuth  cell  of  the  radar  as  the  antenna 
beam  and  platform  both  slowly  change  their  positions.  Since  the 
statistical  system  being  considered  may  be  assumed  to  be  ergodic 
this  time-average  calculation  can  be  replaced  with  a  statistical 
evaluation  of  expectation  over  the  ensemble  of  all  possible 
random  mine  placements. 

There  are  two  common  methods  of  computing  the  ensemble 
expectations:  direct  analytic  evaluation  and  Monte-Carlo  simulation. 
The  direct  analytic  technique  requires  a  probability  density 
description  of  all  the  random  variables  of  the  system  and  becomes 
easily  intractable  for  complicated  systems  having  mutual  dependencies 
between  the  variables.  The  Monte-Carlo  method  requires  a  model  of 
how  the  various  possible  system  configurations  evolve  in  nature 
as  outcomes  of  the  random  process.  The  model  is  simulated 
repeatedly  using  appropriately  distributed  random  number  generators 
for  the  variables  and  expectations  are  obtained  by  averaging  the 
outcomes. 


As  originally  envisioned,  the  expectation  evaluation  was 
to  be  attempted  by  analytic  means  with  the  Monte-Carlo  method  as 
a  backup  in  case  of  failure  In  the  analytic  technique.  As  will  be 
shown,  the  analytical  evaluation  of  the  expectation  can  be  performed 
for  this  problem  using  certain  specified,  but  realistic,  assumptions. 
The  analytic  expression  for  the  expectation  contains  coherent  and 
incoherent  array  factors  which  are  each  sumnations  of  weighted  (by 
the  antenna  pattern)  coherent  and  incoherent  terms  due  to  scattering 
from  the  centroid  positions  of  each  scatterer  in  the  mine  array.  By 
curve  fitting  the  gain  pattern  with  a  decaying  exponential  in 
elevation  and  azimuth  it  is  possible  to  transform  these  two  summa¬ 
tions  into  geometric  series  which  are  then  summed  in  closed  form. 

The  resultant  expression  for  normalized  S/C  displays  the 
inherent  contributions  from  radar  system  parameters  (such  as  beam- 
widths,  frequency,  depression  angle,  and  pulse  width)  in  conjunction 
with  the  randomness  of  mine  placement  and  the  clutter  distribution. 

The  presentation  that  follows  has  been  organized  for 
optimum  useful  information  conveyance  with  the  results  of  this  work 
being  displayed  in  the  next  section  along  with  a  discussion  of 
trends  and  observations.  A  more  detailed  set  of  conclusions  as  well 
as  details  of  the  analysis  that  was  performed  is  left  to  the  final 
section  wherein  a  tutorial  format  is  utilized  to  accommodate  a  wider 
readership.  Some  mathematical  details,  a  special  case,  computer 
programs  and  description  of  notation  are  relegated  to  appendices. 
Finally,  several  references  covering  electromagnetics,  radar 
systems,  radar  cross-section  and  clutter  properties,  probability 
theory  and  scattering  by  random  media  are  included  for  use  by  readers 
and  are  cited  where  appropriate  in  the  text. 


II. 


RESULTS  AND  OBSERVATIONS 


In  order  for  the  results  to  be  meaningful  we  must  first 
consider  the  pertinent  assumptions  and  definitions  Involved  in 
their  derivation.  The  actual  derivation,  however,  will  be  left 
until  Section  IV. 

In  Figure  1  the  geometry  of  the  assumed  model  is  depicted. 

The  radar  system  antenna  illuminates  a  patch  of  the  mine  field  with 
the  centroid  of  the  main  beam  (antenna  boresight)  pointing  towards 
the  origin  of  the  cartesian  coordinate  system.  The  mine  field  is 
assumed  to  extend  over  some  irregular  shaped  region  of  the  x-y 
plane,  n,  covering  many  range-azimuth  cells  of  the  projected  radar 
beam  on  the  surface.  A  range-azimuth  cell  is  depicted  in  part  (b) 
of  Figure  1.  Assuming  that  the  depression  angle  to  the  main  beam 
look  direction  is  0^  and  that  the  range  to  the  surface  centroid 
point  is  Rq,  as  shown,  then  the  surface  projection  of  the  half- power 
contour  of  the  main  beam  will  appear  as  in  Figure  1(b),  where 
0g  is  the  vertical  elevation  (EL)  half-power  (3  dB)  beamwidth  and 
♦g  Is  the  horizontal  azimuth  (AZ)  3  dB  beamwidth,  both  expressed 
in  radians.  The  range-azimuth  cell  is  limited  in  range  by  the 
effective  temporal  pulse-width  of  the  radar,  t.  The  time  gating 
of  the  radar  receiver  allows  range  resolution  of  two  point  targets 

Q 

by  as  little  as  ct/2,  where  c  =  3  x  10  m/sec.  On  the  surface, 
which  is  slanted  with  respect  to  the  look  direction  to  the  radar, 
this  resolution  translates  to  (ct  sec  0q)/2  in  surface  range. 

The  echo  returns  from  both  the  mines  and  clutter  will  be  produced 
from  an  illuminated  surface  patch  with  width  of  several  azimuth 
beamwidths  but  time-windowed  at  one  projected  range  cell.  For 
large  depression  angles  (looking  almost  straight  down)  the  patch 
range  will  be  limited  by  several  elevation  beamwidths  of  the  antenna, 
since  (ct  sec  0q)/2  becomes  Infinite. 


The  mine  placement  Is  statistical.  To  model  this  situation 
on  a  usable  basis  we  assume  that  the  mine  positions  are  randomly 
perturbed  from  some  specific  position  grid,  as  would  be  the  case 
when  mines  are  laid  on  a  systematic  basis  by  moving  from  one 
position  to  another.  The  centroid  grid  will  be  assumed  to  be 
rectangular  with  spacing  of  dx  and  dy.  From  each  of  these 
centroid  positions  a  mine  Is  "thrown"  randomly  with  Its  relative 
X  and  y  coordinates  having  a  zero-mean  gausslan  probability 

density  with  respective  standard  deviations  of  a  and  o  . 

X  y 

The  electromagnetic  echo  from  the  mine  field  is  simply  a 
linear  sum  of  terms  from  each  of  the  random  mines.  This  sum.  If 
containing  a  large  enough  number  of  mines,  can  be  replaced  for 
computational  purposes  with  statistical  expectation  over  the  joint 
probability  density  function  of  the  mine-field  ensemble.  This 
expectation  Is  considered  In  general  and  evaluated  in  closed  form 
In  Section  IV  for  a  special  case  of  practical  significance  where 
the  antenna  gain  pattern  can  be  approximated  In  vertical  and 
horizontal  off-bores Ight  angles  by  an  exponentially  tapered 
function.  A  key  result  from  Section  IV  is  that  the  radar  cross 
section  (RCS)  of  a  mine  field  composed  of  Identical  mines,  each 
having  the  same  independent  RCS^  Is  of  the  form 

RCS„  *  RCS^  .  F  (1) 

where  F  is  an  array  factor  which  Is  a  function  of  the  statistics 
of  the  mine  field  (standard  deviation  of  mine  placements  vs.  wave¬ 
length  of  the  incident  field)  but  Is  Independent  of  the  Individual 
mine  structure  and  RCSM  Is  the  radar  cross  section  of  the  mine  field.  A 
common  assumption  in  taking  expectations  of  random  scatterers  is  that  the 
total  scattered  power  Is  simply  the  sum  of  the  Individual  powers  returned 
from  each  of  the  scatterers,  12;  pp.  77-80] .  This  Is  termed  the  "Incoherent" 
scattering  assumption  and  Is  valid  when  the  standard  deviations  are  much 


larger  than  the  wavelength  of  the  illuminating  field.  On  the 
other  extreme,  when  the  standard  deviations  are  much  smaller  than 
the  wavelength,  the  random  displacements  of  the  mines  from  their 
centroids  do  not  appreciably  perturb  the  phases  of  the  individual 
fields  from  each  of  the  scatterers.  A  good  approximation  is  then 
to  simple  add  the  fields  from  each  of  the  mines  placed  at  its 
respective  centroid.  This  is  termed  the  "coherent"  approximation 
since  the  individual  fields  add  coherently  with  specific  and  known 
relative  phases.  As  is  shown  in  Section  IV  the  array  factor,  F, 
in  equation  (1)  can  be  expressed  as  a  weighted  sum  of  incoherent 
and  coherent  case  array  factors  with  the  relative  proportions 
(weights)  depending  upon  the  "randomness"  of  the  distribution  of 
mines 

2  2 

F  =  e'^  F^  +  (1  -  e"Y  )  Fj  (2a) 

where 

Y  =  Air  ^  cos  8  (2b) 

^0  ° 

for  the  case  of  equal  standard  deviations  a  *  a  =  a  with  A 

X  y  0  0 

the  wavelength  of  the  incident  field  and  the  depression  angle. 
Note  that  as  grows  (increase  randomness)  the  contribution 

of  the  incoherent  term  Fj  increases  relative  to  that  of  the 
coherent  term  F^.  Note  also  that  as  0^  approaches  90®  (look- 
down  incidence)  y  goes  to  zero  and  F  F^.  This  is  true  because 
at  straight-down  incidence  random  perturbations  of  mine  positions 
at  right  angles  to  the  platform  direction  do  not,  under  the  far- 
field  assumption,  change  the  path  lengths  or  phases  of  the  echo 
returns  from  those  of  their  values  for  centroid  located  mines. 

Coincident  with  the  echo  signal  from  the  mines  will  be 
reflected  power  from  the  earth  and  other  objects  colocated  with 
the  mine  field  such  as  trees,  vegetation  and  perhaps  man-made 
structures.  This  unwanted  echo  tends  to  mask  that  of  the  mines  and 


is  thus  termed  "clutter"  due  to  its  cluttering  effect  on  radar 
observation  screens.  The  ratio  of  signal  power  to  clutter  power 
S/C  is  thus  an  important  quantity  for  radar  performance.  As  is 
shown  in  Section  IV  this  power  ratio  depends  upon  the  ratio  of 
mine  field  cross  section  to  that  of  the  clutter  returns 


(3) 


The  effects  of  distance  and  antenna  gain  are  not  explicit  in  this 
ratio  since  these  terms  are  common  to  both  the  mines  and  the 
colocated  clutter  scatterers,  and  thus  cancel. 

It  is  common  to  express  clutter  cross  section  as  a  ratio 
of  cross  section  per  physical  surface  area,  [3].  Using 

this  concept  it  is  shown  in  Section  IV  that  the  signal  to  clutter 
power  ratio  is  given  by 


Z  •  RCsf^  [l  *  -  l)]  <-'> 

Note  that  in  the  limit  of  increasing  randomness  of  mine  placements 
(y  large)  the  S/C  becomes  simply  the  ratio  of  the  individual  mine 
RCS  per  unit  area  (there  is  one  mine  having  RCS^  per  grid  square 
of  area  dx  dy)  to  the  clutter  RCS  per  unit  area.  Th’s  is  the 
incoherent  limit.  For  y  ->  0  (little  randomness)  we  have 


_ 0  c 

dx  dy  .  RCS* 


(5) 


This  is  true  because  the  mines  have  coherent  array  factors  F^ 
while  the  clutter  per  unit  square  has  the  same  incoherent  array 
factor  Fj  as  the  mines  had  under  highly  random  conditions.  Under 
the  condition  of  y  large  the  Fj's  in  the  numerator  and  denominator 
cancelled  leaving  only  the  ratios  of  cross  sections  per  unit  area. 
These  facts  will  be  shown  explicitly  in  Section  IV. 
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The  term  in  curly  brackets  in  (4),  given  by 


will  be  termed  the  signal  to  clutter  "improvement  factor"  since  it 
is  unity  under  the  incoherent  assumption.  Because  of  partially 
coherent  addition  of  random  fields  at  specific  frequencies  for  a 
given  depression  angle  this  improvement  factor  can  become  very 
large--it  can  also  become  less  than  unity  at  certain  frequencies, 
as  will  be  seen.  As  frequency  is  increased  it  ultimately  approaches 
unity  in  all  cases  except  when  the  mine  placement  is  "deterministic" 

(Oo  =  0). 

The  advantage  of  v<riting  the  S/C  in  terms  of  an  improvement 
factor  is  that  all  characteristics  of  the  individual  mines  and 
particular  mine  field  clutter  conditions  are  contained  in  the  cross 
section  ratio  while  the  improvement  factor,  I,  is  independent  of 
the  types  of  mines  and  local  clutter  conditions  but  does  depend 
upon  the  array  size,  randomness,  frequency  and  gain  pattern  of  the 
radar  antenna  over  the  illuminated  area.  The  S/C  is  simply  its 
value  under  "incoherent  conditions"  times  I. 

In  addition  to  the  periodic  frequencies  for  resonances  and 
anti  resonances  of  the  improvement  factor,  to  be  displayed  shortly, 
the  individual  mines  will  have  resonant  behavior  in  RCS^^  over  a 
band  of  frequencies.  This  band  of  frequencies  for  mine  resonances 
will  in  almost  all  foreseeable  cases  be  observed  at  much  higher 
frequencies  than  that  of  the  resonant  band  for  the  grid  improvement 
factor,  I.  At  these  higher  mine  resonant  frequencies  the  improve¬ 
ment  factor  will  be  close  to  unity.  This  topic  is  discussed  more 
extensively  in  the  final  section  (IV.E)  of  this  report  and  is  also 
considered  in  detail  in  a  separate  contract  with  Geo  Electromagnetics. 

To  produce  quantitative  results  several  cases  of  computation 
were  performed  for  the  mine  field  signal  to  clutter  improvement 


factor  given  by  equation  (6).  Two  special  classes  of  geometry  are 
considered.  Class  I  considers  cases  where  the  mine  field  extent  is  large 
compared  to  the  range-azimuth  cell  size  on  the  surface  while  Class  II 
considers  the  opposite  case  of  mine  field  size  much  smaller  than 
the  range-azimuth  cell  size.  In  all  cases  under  both  classes  the 
geometry  was  specialized  to  that  of  Figure  2  with  the  radar  platform 
located  a  surface  distance  of  Rj  from  the  beam  center  ground  point 
while  at  altitude  H^.  A  locally  flat  surface  is  assumed. 

Some  representative  computations  of  improvement  factor  for  the  mine 
field  array  factor  are  shown  in  Figures  3  through  13.  This  set  of  results, 
although  not  exhaustive  by  any  means,  does  indicate  the  general  relation¬ 
ships  between  improvement  factor  as  a  function  of  frequency  and  certain 
key  parameters. 

The  fixed  parameters  in  all  cases  are: 

Rj  =  1000  meters  (ground  distance  from  platform  to 

centroid  of  main  beam  on  the  ground) 

dx  =  dy  =  1.51  meters  (rectangular  grid  spacing  for  centroid 

of  random  mine  positions) 

FREQ.  RANGE  =  10  MHz  to  510  MHz  plotted  in  1  MHz  increments 
(501  points  total) 

In  Figures  3  through  8  Class  I  (mine  field  much  larger  in  extent 
than  the  range-azimuthal  cell  size  of  the  radar)  is  considered  and  in 
Figures  9  through  13  Class  II  (mine  field  much  smaller  than  the  range- 
azimuthal  cell  size  of  the  radar)  is  considered.  The  signal  to  clutter 
improvement  factor  (equation  (6))  is  plotted  versus  frequency  in  Figures 
1  through  8  and  the  array  factor  F  (equation  (2a))  is  plotted  in  Figures 
9  through  13. 
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Figure  4.  Plot  of  I  vs  Frequency  using  same  parameters 


Figure  5.  Plot  of  I  vs.  Frequency  using  same  parameters  as  In  Figure  4, 
except  T  ■  .2  ysec 


Figure  6.  Plot  of  I  vs.  Frequency  using  same  parameters  as  Figure  4, 


Figure  8.  Plot  of  I  vs.  Frequency  using  same  parameters  as  in  Figure  3 
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The  parameters  being  varied,  one  at  a  time,  are  platform 
height  H^,  antenna  azimuth  beamwidth  4>g,  and  mine  placement 
standard  deviation,  a^,  as  well  as  radar  pulse  width,  x.  The 
standard  canonical  or  baseline  whose  spectra  is  displayed  in 
Figure  4  has  parameters 

=  100  meters  (near  grazing  incidence) 

(Jig  =  5  degrees  (narrow  at  the  HF  and  VHF  frequencies 

being  considered) 

=  0.2  meters 

0 

T  =0.1  ysec 

Table  I  summarizes  the  results  of  changing  the  various  parameters 
about  that  of  the  baseline.  Major  observations  are  that: 

(1)  Higher  randomization  (larger  a^)  increases 

attenuation  vs.  frequency  of  resonances.  It  is 
possible  even  with  significant  Oy/dy  ratio  of 
of  0.331  to  observe  a  first  resonance  at 


f 


1 


c 

2  dy 


sec  0^ 


(7) 


(2)  As  indicated  by  (7),  and  Figure  6,  the  resonances 
occur  at  frequencies  where  down  range  centroid 
spacing  dy  is  multiples  of  projected  wavelength 
^0  ®o*  lookdown  angles  produce  higher 

resonances. 


(3)  Increasing  the  number  of  illuminated  mines  in  cross 
range  (wider  azimuth  beamwidth)  increases  resonant 
(coherent)  return  vs.  incoherent  return.  Result  is 
larger  S/C  enhancement  at  resonance  in  proportion 
to  beamwidth  (assuming  narrow  beamwidth). 
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TABLE  1 


FIG. 

NO 


PARAMETER 

CHANGED 


MAJOR  EFFECTS 


3 

a  =  ,01  HI 
o 

a^«  A  over  entire  frea.  range 

produces  highly  periodic  coherent 
resonances. 

4 

Standard 

Case 

Oy  =  0.2m  becomes  wavelength 

comparable  yielding  rapid  decay 
of  coherency  with  frequency. 

5 

T  =  0 , 2  jtisec 

Doubling  the  range  cell  size  also 
doubles  (2N  +  1)  number  of  mines  in 
coherent  function  F^..  Results  in 

higher  "Q"  of  response  function. 

6 

=  lOOOm 

A  lookdown  at  45*  vis-a-vis  near 
grazing  expands  freq.  scale  due  to 
larger  projected  wavelength.  Amplitude 
increase  due  to  larger  range  cell 
projection. 

7 

<p  =  10* 

D 

Wider  azimuth  beamwidth  provides 
larger  surface  area  illumination 
thus  including  more  mines  and 
increasing  amplitude  of  coherent/ 
non-coherent  return. 

8 

a  =  0.5m 

o 

This  causes  a  very  rapid  decay  of 
the  coherency  factor  with  frequency. 

The  primary  resonance  at  100  MHz  is 
greatly  reduced  but  still  highly 
visible. 
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(4)  Increasing  number  of  illuminated  mines  in  down 
range  {longer  pulse  width)  increases  number  of 
interfering  phase  centers  in  the  array  factor. 

This  produces  narrower  bandwidth  of  resonances 
and  si  delobes  in  frequency  and  decreases 
relative  si  delobe  level  thus  enhancing 
observability  of  resonances  in  swept  frequency 
mode. 

In  Figures  9  through  13  the  special  Class  II  of  a  mine¬ 
field  much  smaller  than  the  range-azimuth  cell  of  the  radar  is 
considered.  The  plots  are  of  the  computed  spectra  of  the  array 
factor  rather  than  the  improvement  factor 

F  =  (1  -  e"^^)  Fj  +  e'^^ 

The  standard  case  is  shown  in  Figure  10 

=  100  meters 

=  11  Cross  range  grid  size 
Np  =  9  Down  range  grid  size 

Total  No.  Mines  =  99 

=  0.2  meters 

The  results  are  very  similar  in  behavior  to  that  of  a  large  minefield 
as  per  resonances  and  dependence  on  randomness  of  mine  placement. 

The  S/C  improvement  factor  would  have  an  identical  shape  to  these 
curves  but  would  be  divided  by  Fj  for  the  entire  range-azimuth  cell 
of  tie  radar  as  given  in  equation  (III-l).  The  incoherent  array 
factor  due  to  the  99  mines  is  much  larger  compared  to  the  coherent 
factor  vis-a-vis  the  same  ratio  for  the  previous  set  of  large  mine 
field  results.  The  first  resonance  at  f^  (at  approximately  20  MHz) 
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is  much  less  pronounced  for  even  the  moderately  random  case  of 

=  0.5  m  (Figure  13);  the  lower  sidelobes  of  coherent  interference 
dominate  the  actual  first  resonance. 


Figure  10.  Plot  of  Class  II,  F  vs.  Frequency  using  same  parameters  as  in  Figure  9,  except 
i.e.,  H  =  100  m,  R.  =  1000  m,  =  5",  0„  =  5®,  dx  =  dy  =  1.51  m,  cT^  =  cy^  =  a,,  ' 


Figure  11 


except  with  11  rows  and  9  columns 


Figure  13.  Plot  of  Class  II,  F  vs.  Frequency  using  same  parameters  as  in  Figure  10 
except  =  0.5  m 


III.  METHOD  OF  ANALYSIS 


A.  RADAR  SCATTERING  FUNDAMENTALS 


In  this  section  we  will  briefly  review  the  pertinent 

equations  of  electromagnetic  scattering  that  are  to  be  used  in 

the  mine  field  scattering  case.  Referring  to  Figure  14,  consider 

a  single  radar  "target"  (mine  or  other  object  being  considered) 

which  is  located  at  a  distance  of  R  from  the  radar  which 

transmits  a  total  peak  power  of  W.  watts  into  a  transmit- 

z  2 

receive  antenna  having  gain  function  G(e,^)  =  B  (e,(!))  where 
G  is  maximum  boresight  gain  and  8(0, (|i)  is  normalized 
(maximum  =  unity)  field  directivity  pattern  of  the  antenna.  The 
elevation  and  azimuth  angles  off-boresight  are  e  and  (ji, 
respectively. 

The  incident  power  per  unit  area  at  the  target  location 
will  be  [4] 


i  2 


pi  =  JJL 


4irR^ 


(8) 


where  =  I^Ott  ohms  (free-space  wave  impedance) 

=  incident  electric  field  phasor  (RMS) 

^t’i’t  "  position  in  gain  pattern 

Using  (8)  we  solve  for  the  incident  electric  field  magnitude 


(30  G(0^ 


(9) 


This  can  be  rewritten  in  terms  of  B(0^,4>^)  and  including  the 
relative  phase  of  the  incident  field  phasor  in  terms  of  distance, 
R,  and  wavenumber, 
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i 


Figure  15.  Far-Field  Phase  Path  Length 


The  complex  diagonal  elements  indicate  the  amount  of  self-polarization 
backscattering  while  the  off-diagonals  (reciprocity 

derived  relationship)  indicate  the  relative  cross-polarization 
effects  for  a  particular  aspect  on  the  target.  By  partitioning  the 
incident  electric  field  phasof  into  horizontal  and  vertical  components 


The  RMS  scattered  field  phasor  received  at  the  radar  in  the  "far- 
2D^ 

field"  (R  >  r —  ,  where  D  =  maximum  target  dimension)  is  15] 

*0 
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The  complex  components  of  ^ thus  indicate  relative  amplitudes  and 
phases  of  co-polarization  and  cross- polarization  scattering  for  a 
specified  aspect.  In  general  ^  has  a  very  complicated  behavior  of 
aspect  and  frequency  but  can  be  numerically  computed  for  a  range 
of  simple  object  geometries,  although  often  with  great  computational 
effort. 

To  relate  ^  to  the  usual  radar  cross  section,  RCS^,  of 
the  target  we  note  that  the  power  per  unit  area  returned  to  the 
radar  antenna  in  the  scattered  field  is  14], 

^  irf,2  .  RCS„ 

P®  =  •  P  - 1  (14) 

"o  4.r2 

Using  (8),  (13)  and  (14)  yields 


RCSq  =  4TrR' 


2  l^r 


B.  ARRAY  FACTOR  DERIVATION 

Let  us  now  consider  the  case  of  backscattering  from 
a  random  mine  field  as  previously  described  in  Section  II  and 
illustrated  in  Figures  1  and  2.  We  model  the  statistical  mine 
field  placement  by  independent  random  gauss i an  probability 
densities  with  centroid  positions  at  surface  locations 

=  X-,..  X  +  y„„  y 
mn  mn  mn 

where  the  double  index  (m,n)  signifies  the  grid  number  and  x,  y 
are  unit  vectors  in  the  x  and  y  directions  respectively  (see  Figures 
2  and  15).  This  double  index  is  used  for  the  specific  case  of  a  rectang¬ 
ular  grid  to  be  used  later  in  summing  the  coherent  and  incoherent  array 
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factors,  but  could  just  as  well  be  replaced  by  a  global  grid  number, 
“k."  We  will  use  (m,n)  however,  noting  that  the  grid  arrangement 
has,  per  se,  no  bearing  on  the  expectation  calculation. 

The  S/C  calculation  proceeds  by  defining  this  ratio  in 
terms  of  time-averages  of  mine-array  and  clutter  field  phasor  returns 


S  <eA(t)  <|e/|^ 


The  "A"  subscript  denotes  array  while  "C"  denotes  clutter.  Assuming 

2L^ 

that  the  beam  center  distance  R  »  — r-  where  L=  maximum  dimension 

°  ^o 

of  the  illuminated  patch  and  =  wavelength  of  the  pulse  carrier 
frequency,  then  the  mine  and  clutter  arrays  are  in  the  far-field.  By 
assuming  simple  superposition. 


E  S  .  r  r  E  S 

^  m  n 


(17) 


where  due  to  the  mine  at  the  random  position  J(^  ^  is, 

from  (13), 


r  ^  _  /  f  F  ^  \  ^ 

‘"m,n  ■  ^^m,n  ’  ‘'m,n^  R. 


^^o^m,n 


(18) 


m,n 


and  J  ~  ^  the  same  scattering  matrix  for  each  mine.  All  mines 
m,n 

are  assumed  identical  with  the  same  orientation.  Using  equation  (19), 
the  incident  field  at  the  (m,n)  mine  has  the  form 


eJ_  =  E  ^ 


"■^^o^m.n 


where 


"  '*m,n 

=  [30 


m,n 


(19) 


(20) 


30 


is  the  incident  field  value  at  the  boresight  position  on  the  ground 
(F  =  0)  with  =  Trans  Power,  6^  =  boresight  max  gain,  and  a^ 
incident  unit  polarization  vector,  and 


mn 


mn 


1/2 


(21) 


is  the  normalized  antenna  "field  gain"  at  r  .  The  far-field 

mn 

assumption  allows  replacing  R  „  =  R  in  the  denominators  of  (18) 

mn  0 

and  (19)  while  the  phase  terms  become  6  R  ^  =  e^(R  -  r  •  r  ), 

0  mn  0  0  o  '^n 

as  shown  in  Figure  15.  There  results 


=  (</ 


V- 


Z  Z  B 
m  n 


(22) 


the  signal  power  becomes 


(23) 


Using  (23)  the  equivalent  radar  cross-section  of  the  mine  field  can 
be  easily  found  and  was  given  in  Section  II 


2  <|Ea'|'> 


RCS„  =  4ttR 
n  0 


Thus,  the  average  RCS  of  the  random  ensemble  of  identical  mines 
equals  the  RCS^  of  any  one  mine  in  the  field  times  an  average  array 
factor,  F,  which  is  found  from  (1)  and  (22)  to  be  the  time-average 
of  the  sum  given  by 


F  = 


Z  Z  Z  Z<B„B„.„, 

m  n  m*  n' 


,n'^  y 


(24) 
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The  source  position  direction  unit  vector  is  given  by 


r  =  cos  e  cos  <l>  x  +  cos  6  sin  y  +  sin  z 
0  0  0  0  0  0 


while  the  random  mine  positions  are  given  by 


r  =r  +aUx+aVy 
'^.n  m,n  x  y 


} 


standard 


(25) 


y  "  deviations 

where  the  joint  gaussian  probability  density  is  given  by,  [61, 
P(u.v)= 

where  U  and  V  are  independent  normalized  random  variables. 

Assuming  independence  of  positions  and  system  ergodicity, 

F  =  Z  Z  £[b  ^1 

M  n  ivi*  c^lfnfn 

m  n  tn  n  L 


(26) 


(27) 


(28) 


where 


m  ,n  J 


m  /  m' 
n  /  n’ 


=  cos  0„  cos  (|i„(x„  +  U) 
0  ^11,0  0  0  m  X 


+  cos  0Q  sin  ♦^(y^  +  V) 


(29) 
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Assuming  that  B 


^mn 

i.n  " 


is  essentially  constant  over  several  standard  deviations  in  x  and  y 
from  the  centroid  of  placement  of  the  (m,n)  mine  then 


&  \b  ^1  -  B  ^ 
^  m,n  m,n 


(31a) 


=  B  e 
m,n 


±j23QCOs0^(cos^^  +  sin*^  y^) 


±j26QCOs0^(cos.j>^  +  sin(J>^  V)j 


(31b) 


To  simplify  equation  (31b)  it  is  convenient  to  designate  the 
expectation  in  the  last  line  as  ^  [U  and  V]. 


The  expectation  over  x  and  y  (R.V.s  U  and  V)  becomes  a  product 
of  two  one-dimensional  expectations.  Defining 


“x  '  2  °x 


a.,  =  2  B„cos9  sin*  a 
y  0  0  0  y 


U  and  V 


where 


j  =  S^jon  uj  X  ^l^on  vj 


±ja^  U  -UV2 


(33a) 


(33b) 
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2  2 


Using  J  e"®  ^  cos  bx  dx  "  2a 
0 


.2 
b  =  a 


yields 


Li kewi se 


Thus 


g^on  uj 

<£  jon  vj 
’|lJ  and  vj 


=  e-x/2 


=  e-“y'/2 


^  I^U  and  VJ  =  e 
Substituting  back  into  (28)  yields 

II 


9  9  9  9  9  9 

-23^  cos  0^(a  cos  (f)  +  a  sin  ((>  ) 
0  o'x  ~0  V 


(34) 

(35a) 

(35b) 

(36) 


m  n 


2  -43^  ^  °y^  sin2<J>^) 

be 

m»n 


III  I 


j2S„coseo[coS4.o(Xm  -  x;)  +  -  y^)] 

(37) 


X  La  ^  Kr^ 

~  _i  mn  m  n 

m  n  m  n 

m  ji'  m' 
n  ^  n' 

Defining 

=  43p  005^6^(0^^  cos^4>p  +  sin^(}>p) 

6^  =  23j,coseQ  cos<t>^ 

6y  =  23oCos0q  sw<P^ 

We  can  rewrite  the  total  expectation  in  the  form 

2  r-  2  -v^i  ^  ^ 

B  +  e 

m  ii  '  m  n 


(1  -  e-v^  I  I  B  ^  *  e-v^l  Z  I  B 
mV  ">1  I  V  V 


(38a) 
(38b) 
(38c) 

J(5.x„+6^J.2 


X  m 


in  words. 
Array  factor 


(39) 


incoherent 
array  factor 
(additive  power) 

>^1 


coherent 
array  factor 
(additive  fields) 


34 


(2) 


This  provides  our  equation  (2)  from  Part  II 
F  •  (1  -  F,  t  F^ 

and  mathematically  defines  what  is  meant  by  Fj  and  F^. 

In  this  form  the  relative  amounts  of  power  are  clearly  shown  for  the 

incoherent  and  coherent  factors  of  the  array.  As  is  evident  from 
2 

(39).  if  •y  >>  1  then  the  power  is  essentailly  summed  incoherently 
and  F  -  This  happens  when  the  projected  standard  deviation  of 
the  mine  position  is  much  greater  than  the  wavelength 

(a  ^  cos^4)  +  a  sin^4>  )  cos^0  » 

.  ^  0  0  0,  ^0 

Projected  Variance  of  wavelength 

mine  position  squared 

Square  of  standard  deviation 

2 

Likewise,  if  the  projected  variance  is  <<  x  then  the  power  results 
from  a  coherent  field  summation. 


C.  ARRAY  FACTOR  SUMMATIONS 


The  incoherent  and  coherent  array  factor  summations  can 
be  performed  via  computer  for  any  specified  gain  pattern.  To 
obtain  approximate  results  analytically  we  will  model  (curve 
fit)  the  gain  pattern  in  the  main  beam  to  an  exponential  function 
of  the  form 

I  I 

G(9,(t))  -  G^e  ®  e  ^ 


(41a) 


where 


e  =  elevation  (vertical)  angle  from  bores ight 
=  azimuth  (horizontal)  angle  from  boresight 
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3 


Gg,  <))b  =  3d6  beamwidths 


where 


On 

G{±/,  0)  =  G  (0.±^)  = 


(41b) 


which  requires  that 

a  =  2  In  2  =  1.3863  ^ 

The  absolute  value  signs  in  (41a)  are  put  in  to  make  G(e,(|.)  an  even 

function  of  0  and  By  assuming  that  the  radar  antenna  is  located  in 

the  X  =  0  plane  with  =  -Tt/2,  as  shown  in  Figure  2,  will  allow  an 

easier  evaluation  of  (39)  without  any  real  loss  in  generality  since  our 

coordinate  system  may  be  oriented  as  desired.  The  modified  equations  in 

(38)  become 


6  =  -2B„cos0„ 
y  0  0 

The  projected  gain  pattern  on  the  surface  is  obtained  by  the 
following  approximations  for  a  narrow  beam  (in  azimuth)  and  short 
range  cell 


{38'a) 

(38’ b) 
(38’ c) 


Cross-Range:  l<i>l  -  tan  ^  "  IT"  '*’8 

0  0 


(42a) 


Range-Projection:  -  sin(\-'^ 


-1  y  sin9 
lei  =  sin  — ^ 


for  «  R_ 

sin0p  2  0 


(42b) 
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To  proceed  with  the  evaluation  of  Fj  and  a  grid  geometry  needs 

to  be  specified  for  the  centroids  of  the  statistically  placed  mines. 

We  will  assume  a  rectangular  grid  arrangement  with  x  and  y 

spacings  of  d  and  d  ,  respectively.  The  use  of  a  double  index 
X  y 

(m,n)  anticipated  this  specialization.  Equation  (39)  is  equally 
valid  for  any  specified  centroid  pattern  by  replacing  (m,n)  with 

single  counting  index  k,  with  and  (x^,  y^)  -+  (Xj^,  yj^). 

m  n  k 

The  use  of  a  rectangular  grid  (along  with  the  other  assumptions 
made  up  to  this  point)  will  allow  an  analytic  evaluation  of 
Fj  and  F^.  These  expressions  will  then  yield  the  character! sti c 
type  of  interactions  and  trade-offs  of  parameters  that  are  present 
for  any  type  of  statistical  mine  configuration  being  swept  by  a 
high  resolution  radar.  As  mentioned  previously,  the  case  of  the 
entire  mine  field  occupying  a  small  portion  of  the  range-azimuth 
cell  is  considered  in  Appendix  II. 

Proceeding  with  the  rectangular  grid  assumption: 


m 

-  oc  <  m  <  oc> 

(43a) 

""y 

-N  <  n  <  N 

(43b) 

In  performing  the  summations  of  Fj  and  F^  we  will  limit 
n  by  the  range  cell  size 


crsecGp) 


N  =  Integer 


{43c) 


while  m  is  taken  between  ±<»  because  the  contributions  due  to  non¬ 
existent  mines  beyond  the  actual  mine  field  will  be  weighted  by  the  very 
samll  gain  factors  in  (41a)  so  far  off  boresight.  It  is  not  possible 
to  sum  n  between  ±®  because  in  the  cases  considered  is  small 
and  the  values  of  n  with  |nl  >  N  would  make  a  significant  contribution 
to  the  sum.  From  equations  (21),  (Ala),  (42)  and  (43) 


(44) 


-ad. 


=  e 


^0*^6 


m 


-adySin9p 


substituting  into  (39)  yields 


-ad 

> 

O 


o'^B 


m 


1 1 

E 

n  =  -N 


,  ‘^o®B 


00 


V  hi 


(<15) 


which  is  simply  the  product  of  two  geometric  series  where 


-ad^ 


=  e 


o'^B 


<1  and  Z. 


-adySin9p 


<1 


(45a) 
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1 

i 

i 


i 


These  series  are  easily  summed  using  (I  -1)  in  Appendix  I. 


E 


m  =  -» 


1ml  1  +  Z 


t-tt: 


and 


E 


n  =  -N 


1  +  Z  -  2Z^*^ 
y  y 

1  -  z.'  " 


Using  (38* b),  (38' c),  (43b)  and  (44)  the  coherent  array  factor  in 
(39)  is  expressible  as  a  product  of  two  geometric  series 


fc- 


^  «  ••in 

m  n 


00 

E 

m  =  -00 


-ad. 


o'^'B 


|ml  N 


E 

n  =  -N 


Using  (45' ) 


00 

E 


1ml  N 


'X  /  ^  y 

m  =  -00  n  =  -N 


E. 


n 

Z..  P 


j23„d„cose„ 
where  P  =  e  °  ^  ° 


(46a) 


(46b) 


cose  n  2 


(47a) 


(47b) 
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Using  1-2  in  Appendix  I  we  obtain 


{4Ca) 

=  {1  _  2  .  2(N+1)/2  j-pN+1  ^  (l/P)^'*'^] 

y  y  y 

n  =  -N 

+  [P*^  +  (l/P)’^]}  /{I  +  Zy  -  Zy^  [P  +  (l/P)]} 


=  1  ■-  Zy  -  2{v/^)^''^  COS  [(N+l)2edyCOs0^]  +  cosCNZ^dyCose^] 

1  +  Zy  -  Zjl^  cos(23dyCoseQ) 


While, 

from  (46a) 

with  Zj^-^-Z^, 

yields 

00 

1 

+  v/^ 

E 

zl"'l/2  =  - 

X  1 

-VT 

^  X 

Zx<'' 

m  =  -«> 


(48b) 


The  incoherent  array  factor  can  be  written  in  terms  of  hyperbolic 
trig  functions  as 


ad 


ad  sine 


-(N+ls) 


‘  Vb 


cosh( 


ad  sine 


] 


{45a) 


while  the  coherent  array  factor  may  be  expressed  as 


(4%) 


F 


c 


coth  ( 


ad. 


o'^B 


2 


ad„sin0„  Nad  sInO 

fsinh  (-lor)  -  “P  <■  2R  e„~) 
0  B  0  D 


(N+l)ad  sine 

+  exp  ( - 2re^ - )cos[N2e^dyCos0Q)]}/ 


'o"B 


ad  sine 

{cosh(  2ft  0^  - ) 


-  cosCZe^dyCoseQ)} 


These  expressions  may  be  substituted  into  (39)  to  obtain  an  exact 
formula  for  the  array  factor  which  is  valid  under  the  assumptions 
of  (a)  a  large  mine  field,  much  larger  than  the  antenna  beam- 
range  cell  ground  projection,  and  (b)  almost  constant  gain 

over  several  standard  deviations  of  random  position 

from  each  centroid  in  the  grid.  Certain  realistic  assumptions 
may  now  be  made  to  simplify  these  formulas.  The  first  of  these 
concerns  the  results  of  the  infinite  x-summations  over  m  . 
Assuming  the  gain  factor  exp  |x|^  changes  little  in  a 

single  centroid  step  implies  that  u  =  This  allows  us 

to  truncate  the  coth  term  in  {49b)  at  the  first  term  in  the  power 
series  coth  (u)  =  ^  ^  +  •••  resulting  in 


ad. 


2R 


coth 


d-^B 


2R«<^ 


d'»'B 


ad. 
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(50) 


1 


This  same  result  is  obtained  by  approximating  the  m-series 
in  (45)  and  (47a)  by  continuous  integrals  over  -  «>  <  x  <  «>  . 
The  validity  of  this  integral  approximation  is  also  contingent 
upon  the  small  change  of  the  integrand  with  step  size  in  the 
original  sum. 

Additional  simplifications  result  by  considering  two 
special,  but  realistic,  cases  of  practical  interest.  These 
cases  result  by  considering  the  behavior  of  the  y-projected 
gain  pattern  over  the  range-cell  length  set  by  the  radar  pulse 
width  and  depression  angle.  The  gain  pattern  y-projection 


has  the  form  exp 


(■  ^  ’>') 


where  the  range-cell  has 


ly|  i  ^  sece^. 


Consider  two  cases: 


Case  I 


CTtane, 


<  <  1  implies  that  the  gain  is  approx,  constant 


within  the  range  cell  =  1  and  =  1. 


E 

n  =  -N 


2N  +  1  and  the  incoherent  array  factor 


becomes 


=  (2N  +  1) 


2Ro<^B 


(51a) 
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while  the  Z  sum  in  F  becomes  that  for  a  uniform  array 
y  ^ 


N 


yielding 


F 

c 


|n|/2 


( 


4  R, 


4'b 


ad 


p"  = 

E  - 

n  =  -N 

1 

sin[(2N+l)eQd^cosep] 

/ 

sin(BodyC0s9“)  ■ 

in[(2N+l)3Qdy  cose^] 
sin(3„d  cos0^) 


0  y 


(51b) 


(51c) 


This  case  will,  in  most  applications,  be  the  most  commonly  occurring 
in  a  surface  search  mode  of  operation, with  exceptions  occurring  only 
at  the  closest  ranges  having  nearly  straight  down  depression  angles. 

Case  II 
cetane 

~  a~'b  '  >  >  1  implies  that  the  range-cell  due  to  pulse 

width  is  much  larger  than  the  elevation  beamwidth  projected  gain 
pattern.  In  this  case  we  may  extend  N  to  +  »  in  the  n-series 
for  the  same  reason  that  the  m-series  was  extended  to  +  “  . 

The  end  result  of  this  is  that  exponential  terms  in  (49a)  (49b) 
with  neg.  arguments  containing  the  N  multiplicative  factor 
may  be  considered  to  be  vanishingly  small.  In  addition,  we 
may  approximate  the  hyperbolic  functions  having  small  arguments, 
u  <  <  1,  by  the  first  terms  in  their  power  series: 

•  1  •  *11^ 
coth  u  =  -  ,  sinh  u  =  u  ,  cosh  u  =  1  +  yielding 
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(52a) 


We  are  now  in  a  position  to  complete  the  formulation  of  the 
problem,  which,  restated,  is  to  investigate  optimum  radar  parameters 
for  enhancing  the  signal  to  clutter  power  ratio  in  scanning  the 

mine  field.  This  brings  up  the  topic  of  clutter. 
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D.  SIGNAL  TO  CLUTTER  RATIO 

The  clutter  return  is,  by  definition,  due  to  all 
undesirable  scatterers  being  illuminated  by  the  radar.  Those 
clutter  scatterers  which  are  colocated  with  the  target  in  the  same 
range/cross  range  cell  will  provide  an  echo  signal  which  will  tend  to 
mask  the  target.  In  the  case  being  considered  of  an  airborne 
surface-search  radar,  the  clutter  echoes  will  be  due  to  scattering 
by  the  illuminated  ground  surface  patch  (even  if  the  surface  were 
flat  and  smooth)  plus  scattering  by  any  or  all  individual 
scatterers  such  as  rocks,  boulders,  vegetation  (including  grass 
and  trees),  buildings  and  other  constructed  objects.  The  amounts  of 
scattering  by  the  ground  surface  alone  will  depend  upon  the  grazing 
angle,  frequency  and  degree  of  surface  roughness,  [3]. 

The  analysis  of  total  average  clutter  echo  power  is 
similar  to  that  performed  for  the  case  of  random  mine  placement. 

The  surface  is  modelled  by  an  ensemble  of  totally  random  scatterers 
which  represent  the  undulations  of  the  surface,  as  well  as  grass, 
rocks,  etc.  The  assumption  of  incoherent  power  addition  is  usually 
made.  This  assumption  is  valid  when  the  standard  deviations  of  the 
positions  of  the  scatterers  are  much  greater  than  one  wavelength. 

This  assumption  is  equivalent  to  setting  y  =  ^  in  (39)  which 
makes 

This  assumption  of  incoherence  is  usually  a  valid 
approximation  except  in  special  cases  where  some  periodic  or 
regular  structure  is  present  in  the  surface  clutter  ensemble. 

Examples  would  be  a  plowed  field  or  wire  fence  or  row  of  planted 
trees.  The  coherent  echoes  would  be  enhanced  at  frequencies  where 
phase  addition  or  subtractions  occur— usually  at  periodic  spacings 
which  are  multiples  of  a  projected  half-wavelength. 
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We  win  assume  that  clutter  is  due  to  an  incoherent 
ensemble  of  surface  scatterers  v<hich  are  randomly  distributed 
about  the  range-azimuth  cell  of  the  radar  in  a  homogeneous  manner. 
By  defining  a  continuous  average  cross  section  per  unit  surface 
area  for  the  clutter  scatterers  we  can  use  (1)  in  conjunction 
with  (39)  where  F  =  Fj(y  -»•  “)  to  obtain  an  equivalent  clutter 
cross  section  in  the  range-azimuth  cell. 


RCS^  =  RCS^ 


If 


B  (x,y)  dx  dy 
surface  cell 


(53) 


where  RCSi  =  clutter  cross  section- per  unit  area,  which  is  assumed  to  be 

L  2 

constant,  and  B  (x,y)  =  G(x,y)/G^  =  normalized  antenna  power  gain. 

Using  the  projected  gain  pattern  from  (41a),  (42a)  and  (42b)  and 
assuming  the  range-cell  projection  is  less  than  the  beam  width  elevation 
projection  (see  Figure  lb)  and  that  we  make  negligible  error  in  extending 
the  x  integral  to  infinity 

^  °B  ■'o  ^ 

e  e  dy  dx 


RCSj,  =  RCSJ. 


// 


"  |y| 


r 


2  . 

a  sin0_ 


1 


aCT  tane^ 
%  °B 


(54) 


This  result  can  be  simplified  under  the  special  conditions  of  Cases 
I  and  II  previously  considered, 
c  tan  _ 


Case  I 


<<  1  Gain  constant  in  Range  Cell 


We  replace  e"  =  1  -  u  in  (54)  yielding 

•  ^0  *^B  ^0 

RCS(.  =  RCSJ.  °  ■  -  - - - 


(55) 
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As  mentioned,  this  is  the  case  most  often  encountered. 


CT  tan  e 

Case  II  T-p — r — ^  >>  1  Projected  range-cell  is  large  enough 

that  gain  falls  to  nearly  zero  at  ends  of  range  cell.  Here  we 
ignore  e'“  in  (54),  yielding 

RCS-  =  RCSi  ^  (56) 

^  ^  a  sin 


The  signal  to  clutter  ratio  is  easily  obtained  by  noting 
that  the  echo  returns  from  both  the  mine  field  and  clutter  may  be 
written  in  terms  of  their  respective  cross  sections  by 

G  ^  A  ^  •  (RCS) 

W  =  W.  — ^ - 7 -  Echo  Power  at  Receive  Antenna 

^  ^  (Air)'^  Output  (57) 


where  =  transmitted  power  and  =  wavelength.  In  taking  the 
ratio  (S/C)  all  terms  cancel  except  the  cross  sections  leaving  our 
original  equation  (3)  in  part  II 


S 

C 


RCS,^ 

RCS^ 


(3) 


The  general  formula  for  this  ratio  is  obtained  via  (1),  (39),  (49) 
and  (54)  and  is  computed  using  a  computer  program  described  in 
Appendix  III. 

Rather  than  iterating  the  radar  parameters  blindly  to 
optimize  S/C^in  the  computer  program  for  a  given  case  of  mine 
field  statistics  let  us  step  back  for  a  moment  to  look  at  some 
generic  observations  and  conclusions. 
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E.  GENERAL  CONCLUSIONS 

Our  task  is  to  select  an  optimum  illumination  angle  and  radar 
frequency  so  as  to  maximize  the  signal  echo  from  the  mine  field  while 
minimizing  (or  at  least  not  maximizing)  the  clutter  echo  power. 

Using  (1),  (2),  (3)  and  (54)  the  ratio  we  wish  to  maximize  is 


S 

C 


RCS 


RCSi 


(1  -  e-'f  )  F,  t  e-'<  F(, 

X  - : — X - 


where 


C 


8b  *6 

5  T 

a  sin 


(58) 


1  -  e 


acT  tan  6^ 

^  '*0  8b 


mn 


can  also 


is  the  equivalent  gain-weighted  surface  area  for  incoherent  clutter 
echo.  In  (39)  the  summation  Fj  =  II 

flS  ( 

summation  by  an  integral  as  in  (54). 


be  expressed  approximately  in  terms  of  by  replacing  the 


Fj  =  Ax 


V  // 


B  (x.y)  dx  dy 


(59) 


illuminated 

patch 


where  ax  =  d  and  Ay  =  d  are  the  mine  centroid  spacings 

X  y 

^  weighted  illuminated  area 
"  d  d  ”  patch  area  per  mine 

ad^ 

This  same  result  is  obtained  from  (49a)  when 

^  o^B 

ad  sin  ad  sin  0  -aci  tan  0 

- °«1  and  -(N  +  1/2)  -XL  o_  o 


(60) 


2  Ro 


(61) 


2 

which  is  valid  for  large  N  and  small  changes  of  B  (x.y)  in 
advancing  one  Ax  or  Ay. 
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Dividing  numerator  and  denominator  of  (58)  by  Fj  and  using  (60),  we 
obtain  the  result  discussed  in  Section  II, 


S 

C 


RCS 


'‘x  “c 


X< 


(1  -  e  ^  )  +  e 


!c 


RCS^  1  2  /  F„ 

I'"'''  ' 


(4) 


ratio  of  mine  cross  section  per  unit  surface  area 
to  clutter  cross  section  per  unit  area 


This  is  a  very  useful  result  from  which  to  approach  the  optimization  of 
(S/C).  The  outside  multiplicative  term  is  simply  the  ratio  of  mine  and 
clutter  cross  section  densities.  This  ratio  will  depend  upon  mine 
geometry,  mine  orientation,  polarization,  frequency  and  number  of  mine 
per  unit  area  as  well  as  type  of  clutter  present.  Unless  the  individual 
mines  are  excited  at  a  resonant  frequency  of  the  individual  mine  cross 
section  RCS^  this  ratio  will  probably  be  small.  This  is  one  approach 
which  should  be  investigated  since  the  frequencies  for  resonant  RCS^ 
of  a  typical  mine  will  be  in  the  GHz  regime.  Another  approach  is  to 
optimize  the  apparent  cross  section  improvement  factor  I,  in  the 
curly  brackets.  This  will  be  done  by  first  choosing  a  frequency  so  that 


=  4  cos^  «  1  (62) 

This  requires  a  frequency  low  enough  that  the  projected  wavelength 
is  much  greater  than  the  standard  deviation  of  the  random  mine 
placements  in  the  grid;  A^/cos  9^  »  .  Since  could  be 

several  meters  to  several  dozen  meters  this  could  require 
relatively  low  frequencies  for  normal  radar  operation  (in  the  VHF 
or  even  HF  regions)  or  a  look  down  angle  that  isnearly  vertical 
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W 


(e  -»■  90°).  In  this  latter  case  the  clutter  cross  section  RCSl 
will  generally  rise  due  to  partial  coherence  of  the  surface 
scatterers  ,  [3]. 

2 

Assuming  y  1  then  the  frequency  and  aspect 
angle  must  be  chosen  to  increase  »  Fj  .  This  can  be  done 
by  selecting 

B^d^  cos  =  nrr  (n  =  integer)  (63) 


so  there  is  phase  coherence  between  the  echoes  of  the  centroid 
positions.  The  simultaneous  solutions  of  (62)  and  (63)  will 
only  be  possible  if 


-  << 


1 

2n7r 


(64) 


1 

ideally--say  -^  <  t  instance.  This  restriction  states 
■'  dy  '  471 

physically  that  the  randomization  (o^)  of  the  grid  should  be 
much  smaller  than  the  average  grid  spacing.  This  ratio  indicates 
how  random  the  mine  placement  is.  If  the  grid  is  too  random  then 
F  =  Fj  and  we  are  stuck  with  trying  to  resonate  the  echoes  from 
individual  mine  scatterers. 
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APPENDIX  I.  SUMMATION  FORMULAE 

The  summations  to  be  performed  are  of  the  forms 


k,  (kp+l) 

Z  ^  -  2 

|kJ+l  (k.  +  l) 
1  +  z  -  z  ^  -z 

Ikol  Ikiln 

z  -  z 


k2  >  kj  ^  0 

k2  2-  0  >  kj  (I-l) 
0  >  k^  >  kj 


k^  >  kj  >  0 

^2,  —  ^  ^ 

0  >  k^  >  kj 


(1-2) 


1 

I 

I 

1 


APPEND  X  II.  SCATTEIIING  FACTORS  FOR  SMALL  MINE  FIELDS 

For  the  case  of  a  mine  field  which  occupies  only  a  small 
portion  of  the  surface  projected  area  of  the  radar  range-azimuth 
cell  certain  modifications  need  to  be  made  in  evaluating  Fj 
and  F.  for  insertion  into  (39)  and  (58).  This  special  case  could 
occur  when  a  wide  beamwidth,  long  ranges,  or  small  mine  field 
areas  are  encountered. 

The  assumption  here  is  that  the  gain  value  on  all 

2  2 

illuminated  mines  is  essentially  constant  with  B  =  B  <  1 

m,n  0 

depending  upon  the  mine  field  position  in  the  beam.  From  (39) 
we  then  have 

(II-l) 


(II-2) 


in-3) 

The  geometric  series  in  (II-2)  provides  a  result  which  is  identical 
to  that  of  the  array  factor  for  a  uniform  linear  array  with  element 
spacing  of  2  dy.  (See  (47b)  and  (51b)). 


F^  =  YY.  B  (2M  +  1)(2M  +  1) 

i  nin  0  N - ^ 


m  n 

From  (38' ) ,  (39)  and  (43) 


fe¬ 


ll 


m  n 


total  §  mines 
2 


B  e 
mn 


-jnBBgdyCose^ 
N 


=  (2M  +  1)^  B.J- 


1  P" 

n=-N 


where 


p  =  e 


j2BodyCOS0^ 
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(II- 


F(,  =  (2M  +  1)^ 


sin[(2N+  l)B^dyCOse^] 
rintB^dy  cosGq] 


When  the  dimension  (2N  +  1)  dy  »  this  coherent  array  factor  will 
be  a  very  rapidly  changing  function  of  6^.  This  array  factor  will 
have  a  lobe  structure  with 


Fp  =  (2M  +  1)^  (2N  +  1)^  B  ^  (II- 

^max 

when 


and  2N  -  1  secondary  maxima  about  13.5  dB  below  F^. 

^max 

approximately  equispaced  in  9^  between  global  maxima.  There 
are  also  2N  nulls  located  between  each  Fj.  angle  -  again  with 
approximately  equal  spacing. 

The  lobe  structure  in  F^  will  be  very  fine  for  a  mine 
field  dimension  (2N  +  1)  »  X^  with  angles  between  nulls  of 


Ae 


null 


2N  d. 


radians 


(II- 


For  Example;  f  =  10  GHz  and  2N  d„  =  20  meters  (range  dimension  of 


y  mine  field) 

yields  Ae  , , »  =  .039  rad  s  2.2“  (spacing  between  nulls] 


APPENDIX  III.  COMPUTER  PROGRAMS 

Case  I.  Programming  (S/C)  Using  Equation  (58) 

Using  (59)  which  assumes  that 


1 

O 


X  y 


// 

illuminated 

patch 


B  (x,y) 


(III-l) 


acTtane, 


The  replacement  of  the assumes  that  B^(x,y)  has  small  changes 

from  one  mine  centroid  to  the  next  and  that  the  mine  field  extends 
far  beyond  the  beamwidth  in  the  transverse  direction  yielding 
-CO  <  m  <  «.  Using  this  same  assumption  the  coherent  array  factor 
summation  in  m  can  be  converted  to  an  integral.  The  result  of  this 
is  obtained  by  the  approximation  in  (51c) 

.2 


sin  (2N  +  1)  Bpdy  cose^ 

[ 

sin  (e^dy  cos0p) 

(I1I-2) 


where  (2N  +  1)  =  Integer  Part  of 


CT  sec0j 

2  d.. 


=  No.  of  mines  within  the  pulse  range 
resolution  in  the  y-di recti on 
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This  approx,  also  assumes  that  the  gain  does  not  change  significantly 
as  a  function  of  range  within  the  range-cell  projection: 


cxtane^ 

~4R^e^  <  <  1  Assuming  low  grazing  angles  (III-3) 


The  equation  to  be  programmed  is  (4)  in  Section  III-E 


RCS„  1 

1  ' 

> 

"  d"d  j 

''x^y  C  1 

mine  to  clutter 
cross  section  density  ratio 


where  y  =  ^^o'^y 


(III-5) 


The  resultant  computer  code  as  programmed  in  BASIC  for  use  on  an 
HP-85  microcomputer  is  shown  in  Figure  III-l.  Input  parameters  as 
requested  by  the  program  are: 


Platform  Height  (meters) 

Ground  distance  to  Main  Beam  Control 

Vertical  and  Horizontal  Beamwidths 
in  Degrees 

AX  (cross  range)  and  aY  (down  range) 
grid  centroid  spacings  in  meter 
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Fioure  III-l 
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Large  Mine  Field  Case 
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Case  II.  Special  Case  Prog^ram  -  Uniform  Illumination 


The  total  Array  Factor,  fran  (39)  and  (38'a),  is 
2.  2 

F  =  (1  -  e  F,  t  e  -If  F^ 


p  ?  2  2 

where  y  =  ^3^0^  cos  0^ 


is  the  coherency  factor  with  o^  the  std.  deviation 


Assuming  B 


m 


- 

.0  “  I  ^  ' 


li 


1  all  m,n 


for  the  case  of  uniform  boresight  illumination  we  have  from 
(II-l) 


Fj  =  (2M+1)-(2N+1)  total  No.  of  Mines 


where  2f1+l 
and  (II-4) 


cross  range  columns,  2N+1  =  down  range  rows 


=  (2ri+l)^ 


sin{(2N+l)8Qdy  cose^) 
sin  (8_d..  cose.) 


The  equivalent  RCS  of  the  array  is  from  (1) 


RCS„  =  RCS„  F 
M  0 


(III-ll) 


assuming  uniform  mine  shapes,  composition  &  orientations 


The  programming  of  F  from  (III-7)  for  this  special  case  of  uniform 
illumination  (Appendix  II)  results  in  the  BASIC  code  shown  in 
Figure  III-2.  This  code  is  a  modified  version  of  that  in  Figure 
III-l  with  inputs; 


Platform  Height  and  Ground  Distance 

Cross  Range  Number  of  Columns  in 
the  Basic  Mine  Grid 

Down  Range  Number  of  Rows 

AX,  AY  Grid  Centroid  Spacings 

Oy  (Standard  Deviation) 

Start  and  Stop  Frequencies  (GHz) 

No.  of  Frequency  Points 


Output: 

_  2  _  2 

F3(I)  =  (1  -  e"^  )  Fj  +  e~^  %  ”  ^  "improvement  factor" 

F5  <  f  <  F6 


for 
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Picture  III-2  BASIC  Computer  Program 
Small  Mine  Pi  el  cl  Case 
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